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f— ^ ■ We study the discovery potential of the Minimal Universal Extra Dimension model 

(MUED) and improve it utilizing the multijet + lepton mode at the LHC. Since the 
MUED has a nearly degenerate spectrum, most events only have soft jets and small 
i?™**. The signature is challenging to search. We apply Mt2 for the event selection 
K> , and set the invisible particle mass of Mt2 (test mass) to zero. The test mass is much 

^ I smaller than the invisible particle mass of MUED. In that case, Mt2 of the signal can be 

large depending on Up-Stream Radiations (USR) which includes initial state radiations 
(ISR). On the other hand, Mt2 of the background is mainly below the top quark mass. 
Hence, the signal is extracted from the background in the high Mt2 region. Since we 
use the leading jets for Mt2, there is a combinatorics effect. We found the effect also 
enhances the signal to background ratio for high Mt2- We perform a detailed simulation 
with the Matrix Element correction to the QCD radiations. The discovery potential of 
the MUED is improved by the Mt2 cut, and especially, the improvement is significant 
for the most degenerate parameter we consider, AR = 10. 
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1 Introduction 

In the last decade, extra dimensional models were studied as possible extensions of the 
Standard Model (SM). The Universal Extra Dimensions model (UED) was developed by 
Appelquist, Cheng, and Dobrescu [T] (see [2] for review). While other types of extra dimension 
like the ADD model [3] and the RS model [1] are meant to solve the gauge hierarchy problem, 
the UED is motivated by the Dark Matter problem [5]. One of the most attractive ways to 
explain Dark Matter is a new weakly interacting particle (WIMP), and the UED contains 
WIMP. 

In the UED, all the SM fields propagate in the flat extra dimensions. The fields are 
expanded in discrete Kaluza-Klein (KK) modes called KK particles according to the KK 
number, the extra dimensional momentum, in the 4D effective theory. In the 5D UED, the 
extra dimension is compactified on an S^ jZ^ orbifold to obtain the SM chiral fermions in the 
zero mode. This orbifold violates momentum conservation in the extra dimension but KK 
parity remains unbroken. KK number odd (even) modes have odd (even) KK parity, and the 
SM particles have the even KK parity. The lightest KK odd particle (LKP) is stable because 
this cannot decay into the SM particles due to the KK parity. 

We consider the Minimal Universal Extra Dimension model (MUED). The mass spectrum 
of each KK level at the tree-level is highly degenerate in mass: mass of the nth KK level 
is approximately n/i?, where R is the radius of the compactified extra dimension. The 
degeneracy is a little relaxed due to radiative corrections at the one-loop level [6l [7] . The 
contribution of the radiative corrections is large when KR is large, where A is a cutoff scale 
of the MUED. However, KR cannot be larger than KR ~ 40, because the running gauge 
coupling of f/(l)y increases as power law beyond the MUED scale \jR and blows up at 
~ AQ/R. So the mass spectrum is still nearly degenerate. 

For most choices of parameters, LKP is the first KK photon 7*^^^ which is a good Dark 
Matter candidate. In calculating the relic abundance of LKP 7*^^^ in the nearly degenerate 
spectrum, the co-annihilation effect takes an important role [HI [H [9] . The second KK particles 
also enter in the computation at the one loop level [lOl [TTl [12] . The mass scale of LKP Dark 
Matter consistent with the cosmological observations is 1/ R ~ 1.5 TeV [T2] . 

The collider signature of the model with the nearly degenerate spectrum is more difficult 
to find because produced visible particles tend to be soft [13]. The well-studied new physics 
model is Supersymmetry (SUSY) with large mass splittings, and the decay of the colored 
SUSY particles produces hard jets and large E'™*^. However, in the MUED, soft jets and 
small i?^**'^ are mostly generated due to the spectrum. 

The discovery studies at Hadron Colliders have been carried out since the MUED was 



proposed [lU [13 [ISl [ITJ dS] . Previous studies were mainly based on the multilepton channels 
[T4l [T6l \T7\ [T8] and the most promising one is 41 + Eip^^^ in which the background level is 
quite low [HI [161 [I3- But the signal events remain very small in this analysis because this 
channel is accessible to very limited production processes of MUED. The analysis including 
background systematic uncertainties was studied by the CMS collaboration [16], and the 
discovery reach is 1/ R ~ 600 GeV with 1 fb~^ at ^/i = 14 TeV. Since the LKP dark matter 
scenario favors the very high mass scale miKP — ^/R ~ 1-5 TeV [12], it is challenging to 
look for the signature in Al + EJp^'^'^. In order to check the scenario at the LHC, the sensitivity 
should be improved, and therefore, we need an alternative way to search for the MUED. 

Multijet channels without requiring multilepton have a statistical advantage because most 
MUED events have multiple jets, even though the analysis based on hard jets and large £^™^* 
cannot deal with the signal of the nearly degenerate model due to the softness of jets and 
small £'™'^*. The problem of the ordinary multijet +E^^^^ analysis is that the signal is buried 
in the SM background. This is because the top quark pair production ti generates missing 
particles, neutrinos, and hard jets with a large enough cross section. If there exists a method 
that can extract the signal from the background based on jets, the discovery potential of the 
MUED could be improved. 

We tackle the search in the multijet channel by using a kinematic variable Mt2 [El [20] . 
sometimes called "Stransverse Mass". Mt2 was originally proposed to measure masses of 
pair-produced particles in the situation with two invisible particles. When the true mass of 
invisible particle is given, Mt2 is bounded by the mass of the produced particles. It was 
proposed that Mt2 can be used not only as a mass measurement variable but for the event 
selection [2T], [22], and this has been applied for the SUSY search by the ATLAS [22] and 
CMS [21] Collaborations. 

In this paper, we point out that Mt2 is effective for the search of the nearly degenerate 
model like the MUED. To use Mt2 as an event selection, we need to set a test mass for the 
invisible particle, and it is set to zero. The set test mass is wrong for the MUED events 
because it is much smaller than the mass of LKP. This leads to the Mt2 dependence on Up- 
stream Radiations (USR). USR is defined as visible particles which contribute to the recoil 
momentum of the subsystem of the pair-produced particles, and they are mainly initial state 
radiations (ISR). Mt2 of the signal can be large depending on USR, although, without USR, 
Mt2 is small in the nearly degenerate spectrum. On the other hand, the test mass is correct 
for the mass of the SM invisible particle, neutrino. Then, Mt2 of the SM background does 
not depend on USR. As shown in Refs. [191 [20], it is mainly below the mass of the heaviest 
particle in the SM, the top quark, Mf^ < mtop- Therefore, an excess in the high Mt2 region 
beyond mtop can be seen as the signal of a nearly degenerate model, and then, Mt2 is effective 



to search for the MUED. 

In the analysis of this paper, leading two jets in pt are used to calculate Mt2- They do 
not always correspond to jets we want, that is, we have combinatoric issues when choosing 
jets for defining Mt2- Combinatorics smears the Mt2 distribution, and the smearing effect is 
different in each process. We found combinatorics makes Mt2 of the signal larger while M7-2 
of the background does not increase as much as that of the signal. Therefore, the smearing 
effect of combinatorics enhances the signal excess in the high M7-2 region. 

We apply Mt2 to the discovery study of the MUED, and we require at least one lepton 
in addition to multijet to avoid the QCD background. Since the ISR takes an important role 
in this method, we perform the event generation with the Matrix Element correction which 
evaluates the hard ISR appropriately. This way, we show that the Mt2 analysis improves the 
discovery potential compared to the 41 + EJ^^^'^ analysis. 

The paper is organized as follows: in section 2 we briefly review the MUED and describe 
its LHC signature, the relic abundance of the LKP, and the experimental constraints. In 
section 3, we discuss Mt2 for the event selection in searching for the MUED signal. The 
simulation setup is presented in section 4. We present our analysis result and show that the 
discovery potential of the MUED is improved in section 5. Section 6 contains discussion and 
conclusion. 

2 The Minimal Universal Extra Dimension model 

2.1 Setup 

In the case of the 5D UED, there is a compactified fiat extra dimension in which all the SM 
fields universally propagate in addition to the 4D Minkowski space-time. Fields are expanded 
in the KK modes (KK particles) in the 4D effective theory, and each zero mode corresponds 
to the SM particle. For example, the 5D real scalar field is decomposed in an infinite number 
of the KK modes after integration of the compactified extra dimension y. 
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and m^ = m|^/ + {n/Ry. R denotes the radius of the extra dimension, and tusm denotes a 
SM particle mass. The fifth dimensional momentum is the mass in the 4D effective theory, 
and this is much greater than nisM, because 1/R ~ O(TeV). Therefore, we can neglect msM'- 
rrin — n/ R, which means the mass spectrum of each KK level is highly degenerate. 

Since the simple compactified extra dimension S^ gives vector-like fermions, an orbifold 
compactified extra dimension S^ jZ^ with an identification of y ^ —y is considered in order 
to obtain chiral fermions in the zero mode. The orbifold compactification results in another 
significant characteristic, the KK parity. KK number is conserved by virtue of the fifth 
dimensional momentum conservation on S^ compactification, but this is broken down to the 
KK parity by the orbifold compactification. The KK parity refiects "evenness" and "oddness" 
of the KK number. All the SM particles have the even KK parity. The lightest particle with 
the odd KK parity, called the lightest Kaluza-Klein particle (LKP), is stable since it cannot 
decay into lighter SM particles due to its oddness. The stable LKP, typically the first KK 
photon 7*^^^ can be a weakly interacting massive particle (WIMP) and therefore a good Dark 
Matter candidate. 

To discuss collider phenomenology, we have to determine the mass spectrum. In this 
paper, we discuss the Minimal Universal Extra Dimenison model (MUED). The MUED is a 
minimal extension of the 4D SM Lagrangian to the 5D UED. At the cutoff scale A it contains 
only SM fields and no other terms, especially no localized terms at two fixed points y = 0, vri? 
led by orbifold compactification. The model parameters of MUED are only three: 5D radius 
i?, cutoff scale of MUED A, and the SM Higgs mass rtih. 

1.1 Mass spectrum 

Radiative corrections to masses of the KK modes at the one-loop level were studied in Refs. O 
[7]. This correction enlarges mass splitting for each KK level away from the highly degenerate 
mass spectrum. The corrected masses are: 

77/ 

"^x(") = ^ + "^x(o) + "^"^xC") (Fermion), (5) 

where ?Tiy(o) is a SM particle (zero mode) mass. Since XjR is taken to be larger than 400 GeV 
as mentioned in Sec. 12.41 Tn^(o) is much smaller than \jR. The neutral gauge bosons of U[X)y 



and SU{2)l are mixed up in the SM, but mass eigenstates of the KK neutral gauge bosons, 
7*^"^ and Z'^'^\ are nearly U{1)y and SU{2)l gauge eigenstates, -B'-"-* and W^^^\ respectively 
because the diagonal components of mass matrix dominates as 

where g' is the gauge coupling of t/(l)y, g is that of SU{2)l, and v = 246 GeV is the vacuum 
expectation value of the Higgs field. The radiative corrections to gauge boson masses are 
given by 
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where C(3) = 1.20205... and gs is the gauge coupling of SU{?>)c- The second terms in the 
corrections are dominant, so mjy{„) is lifted, and m^{n) is slightly lowered. 

The mixings of the KK quarks (KK leptons) are also negligible, and they become f/(l)y 
and SU{2)l gauge eigenstates. Neglecting msu, radiative corrections to the KK quarks and 
KK leptons are given by 
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where Q(") and L^") denote the SU{2)l doublet, and m^"), S''\ and e^") denote the SU{2)l 
singlet. For the KK top quark we should consider the correction from its Yukawa coupling, 
but the production cross section is small. We do not consider the processes of KK top quark. 
Most KK particles receive positive mass corrections. The heaviest particle in each level 
is g^^'^ for the largest correction, and the lightest particle in each level is typically 7^"^ Then 



the LKP is 7^^^ with the mass m (i) = 1/ R. If the Higgs boson of the SM is as heavy as 
rrih ^ 240 GeV, the first KK charged Higgs h^'^^' can be LKP due to the negative mass 
correction of the Higgs four point couphng. But, of course, this cannot be the Dark Matter. 
In this paper, we fix the Higgs mass at m^, = 120 GeV, and we keep 7'-^^ as LKP so that it 
is the Dark Matter candidate. 

The corrections are basically proportional to InAi?, so the degeneracy is crucial for the 
smaller KR. The cutoff scale of the UED was discussed in [23], and the appropriate cutoff 
scale should be several dozen \/R for a given R. As the energy scale grows, more KK particles 
appear, and the logarithmic running of the gauge coupling changes into power law running 
above the MUED scale 1/ R. The f/(l)y gauge coupling blows up (Landau pole) at the energy 
scale ~ 40/i?, so we should set the cutoff scale below the Landau pole. The very small KR is 
also not appropriate because we should consider the higher dimensional operators, and the 
MUED framework is not a good effective theory any more. In our analysis, we considered 
10 < Ai? < 40. A benchmark point of MUED is chosen as l/R = 800 GeV, AR = 20, and 
table d] shows its mass spectrum. 
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Table 1: Mass spectrum of first KK level for a benchmark point {l/R, AR) = (800, 20) 



2.3 Production and decay at the LHC 

At the LHC, the first KK particles of the odd KK parity are pair-produced, and they even- 
tually decay into the LKP. The dominant production processes are KK gluon-|-KK quark 
{g^^^ + Q^^^q^^^) and KK quark+KK quark (Q^^V?^^^ + Q^^V?^^^)- The cross sections of the 
colored particles are shown in [13]. For our benchmark point, a{g^^^ +Q^^y q^^^) = 12.2 pb and 
a(Q(iV?^^^ + Q^^V?^^^) = 7.4 pb at v^ = 14 TeV. The g'^^^ decays into Q^^^Q and g^^^g with 
branching ratios, BR{g — )■ Q^^^Q) ~ 40% and BR{g — )■ q^^^q) ~ 60%, respectively. The ratio 
of inclusive KK quark productions is roughly Q^^^Q^^^ : q^^'^q^^'^ : Q^^'^q^^^ = 1:1:2. Because 
q^^'^ only has the U{1)y gauge interaction, it directly decays into 7^^^g. The hard jets mainly 
come from this decay. The branching ratios of Q^^^ are typically BIl{Q^^^ — > QW^^^^) ~ 65%, 
BR(g(^) -^ QZ^^^) ~ 32%, and BR(Q(i) -^ Q-^^^^) ~ 3%. Once W^^^ and Z^^^ appear from 
Q*-^-*, they cannot decay hadronically for kinematical reasons. They democratically decay 
into all lepton flavors: VT^^^) — )• ^^^^lu and Z^^^ — )■ 'j^^^uu or 'y^^H^l" through l'^^^ or z/*^^\ 

This collider signature has been studied in clean channels of multilepton such as 4/ + 
j^miss p^ [T^ [T7] . dilepton, and trilepton [T71 HB]. The missing transverse energy EJp'^'^^ is 



a magnitude of missing transverse momentum p™''* which is measured by the negative sum 
of transverse momenta of visible particles , p™'^^ = — ^p^^ — ^Pt^- The leptons arise 
only from the KK gauge boson W^^"^ or Z^^'^ production. The 4/ + E™'^^ channel has been 
the most promising one because the background is extremely small, but the fraction of the 
MUED events going to this channel is about 1%: from the Q^^^Q^^^ production, each Q^^^ 
should decay as Q^^^ — )■ QZ^^^ — )■ Ql'^l'-y^^^ with the branching ratio of 16%. 

Multijet channels without requiring multileptons are statistically advantageous, so we use 
the multijet + lepton channel. This is accessible to about 65% of the MUED total production. 
The requirement of one lepton is only to avoid the QCD background. 

However, we face the difficulty of relativity small p^^ and £'™'^* due to the small mass 
splitting between produced colored particles and LKP. Then, the ordinary multijet + EJp^'^'^ 
analysis optimized for the typical SUSY expecting hard jets and large E™^'^ is not effective 
for the search for the nearly degenerate model: MUED. It is important to study a way to 
squarely address the nearly degenerate model using multijet because the discovery potential 
of the MUED could be improved for the statistical advantage. In order to search the signal 
of the nearly degenerate model in the channel we choose, we apply Mt2 to an event selection 
described in Sec. 3. 

2.4 Relic abundance of LKP and experimental constraints 

The LKP Dark Matter relic abundance has been studied in Refs. [HI [H [9l [TOl [H], [12] . Inclusion 
of co-annihilation processes is essential for the correct estimation of the LKP mass. In the co- 
annihilation calculation, the effect of s-channel resonances of the second KK particles [TUl [H] 
and the effect of the final states with the second KK particles [12] must be considered. The 
MUED mass scale 1/R consistent with cosmological observations is 1/ R ~ 1.5 TeV including 
these effects [12]. In this case, it is challenging to discover the MUED both at the first 
run (-y/i = 7 TeV) and at the second run (a/s = 14 TeV) with a low luminosity 0(1) fb~^. 
Hence, it is important to enhance the discovery potential of the MUED by developing a new 
technique and/or a new channel to check the LKP Dark Matter scenario at the LHC 

The electroweak precision test suggests the energy scale of the extra dimension 1/ R should 
be greater then 550 GeV for mh = 100 GeV and 350 GeV for rrih = 400 GeV at the 95% 
C.L. [261 [23 • Independent of the Higgs mass, the observed branching ratio of B^ — )■ ^^7 
constrains 1/R > 600GeV at the 95% C.L. [281 [29l [30] . 

We have obtained a collider bound on the MUED from the ATLAS SUSY search in 
multijet +E^'^^^ with 35 pb~^ [23] and 1 fb~^ [31] data and in multijet + one lepton +E^^^'^ 
with 1 fb^^ data [32]. It gives an upper limit on cross sections after several cuts. We have 



checked which parameters of MUED are excluded using Monte Carlo samples generated by 
Pythia[33] 0, and we found that l/R < 600 GeV with 10 < A < 40 is excluded by the 
multijet +EJ^^'' analysis [3lj at the 95% C.L. 

Then, we focus on 400 GeV < l/R < 1600 GeV due to the LKP abundance, the branching 
ratio of B^ — )■ ^^7, the elctroweak precision test, and the LHC constraint. 

3 Method for Searching for the MUED 

3.1 Mt2 

Our idea is to apply Mt2 to the event selection when searching for the MUED. The effec- 
tiveness of Mt2 to search for SUSY were discussed in [21], [22] , and Mt2 was already applied 
in the search for SUSY in multijet + EJ^'"' by the ATLAS collaboration [23]. 

We briefly review the definition of Mrp2- Mt2, an extension of transverse mass M^, 
was originally proposed as a mass measurement variable in the situation with two invisible 
particles [191 [20]. In each event, we only know the total missing transverse momentum, p™^*, 
but each transverse momentum of the invisible particle cannot be measured. The definition 
of Mt2 is : 
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where My is defined by 

Mrp — MT{mms^i),minv{i),PT ,Pt 
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The transverse energy Erp is given by 



Et = ym2+|p^|2. (11) 

In calculating Mt2, we first construct transverse mass M-f,*~ ' and take the maximum of 
them for one partition of p™^ and p^'' satisfying p^'' + p^'' = p™*^. Then, all the 
possible partitions are considered, and the minimum value among them is taken. 



"'^The event generation of MUED that was used for this bound is described in Sec. 4.1. 
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Figure 1: Schematic picture of typical MUED process, g^^'^q^^'> production. When two jets 



from q^ ' — )■ j^ 'q and used as visible particles to construct Mj'2, a jet from the decay of 
^(1) ^ qWq and ISR become USR. 



Let us consider the simple case where the same parent particles are produced and each 
of them directly decay to a visible particle and an invisible particle. If the invisible particle 
mass ruinv is known, Mt is bounded by the parent particle mass, Mt < rriparent in the correct 
partition. Then, as seen from the definition, Mt2 is also bounded by the parent particle 
mass, Mt2 < i^parent- Wc cau also interprete M^2^ as the invariant mass. 

In practice, rriinv is not known. In calculating Mt2, we need to set a test mass for the 
invisible particle. Many attempts were made to simultaneously determine the masses of the 
parent and invisible particles [511 133 ESI EZl EH EH]- One of them [32] utilizes the effect of 
Up-stream Radiations (USR). USR is defined as visible particles which are emitted before 
parent particles of our interest are produced. The transverse momentum of USR, P-^, is given 

by, 



vis(l) , vis(2) , rni 

Pr + Pt + Pt 



(USR). 



(12) 



Pji is a measure of the recoil of the parent particles. The source of USR is mainly initial 
state radiations (ISR). The decay products can be USR if the decay is before the production 
of the parent particles. Fig. 1 illustrates USR in a case of g^^^'q^^'' production. When we are 
interested in q^^\ a quark emitted from the decay of g^^^ and ISR are considered as USR. 

Of course, M^2^ has different behaviors depending on whether the test mass is correct. 
When we set a correct test mass, M^2^ corresponds to the parent particle mass independent 
of USR. But, when we set a wrong test mass, M^2^ varies with USR. This is because Mt 
varies with USR and is no longer bounded by the parent particle mass. This property can 
be used in the search for MUED. 



3.2 Mt2 for the event selection 

The features of Mrp2 for the purpose of event selection were studied in pTl [22] • We use the 
Mt2 cut as an event selection setting the test mass to zero. The set test mass is a correct 
value for SM because the only invisible particles of the SM are neutrinos. In this case, Mt2 
of the most background events, especially ti events, is lower than the top quark mass ?7itop- 
This is because we can measure the parent particle mass with the correct test mass, and the 
top quark is the heaviest parent particle in the SM. Also, it was found that events without 
missing momentum or with fake missing momentum which is parallel to a mismeasured jet 
have very small values of Mt2 [21] • If a significant excess of M'r2 above rritop is observed, it 
should be the new physics signal. 

On the other hand, Mt2 of new physics is not bounded by the parent particle mass because 
the test mass is wrong for the Dark Matter candidate of new physics. The upper bound of 
Mt2 is a mass combination of the parent particle and the invisible particle in the absence of 
USR, 

2 2 

777 — 777 

"^^parent 

In this case, the signal is extracted from the background for models with a large mass splitting, 
such as SUSY, but not for models with a degenerate mass spectrum, such as yUo < ''^^top, 
because the signal is buried in the background. 

However, considering the recoil momentum of parent particles by USR, Mt2 is still a 
useful variable for the event selection in searching for the nearly degenerate model: MUED. 
M^2^ varies with USR and can exceed /io due to the wrong test mass. For example, when 
parent particles of same mass are produced and directly decay to invisible particles emitting 
visible particles (g(^)g(i) -^ qqi^^h^^^), M^^ [3H1 EH] is 

M^r = VMi¥FTi^^^I(^ > /^o, (14) 



where 

MPt) - /^oh/i+lTe:^^) -^T^\ (15) 

'^'ii'parent 




-^ Ho -^ for Pt > niparent 



where Pt is the magnitude of the momentum of USR. There is a rich source of USR because 
processes of heavy particles tend to come along with hard QCD radiations including ISR. 
Hard ISR gives large recoil of produced particles, and M^2^ can have a large value depending 
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on USR. Note that the background events do not have Mj'2 dependence on USR because the 
test mass is correct, so most events are kept lower than mtop. 

When analyzing events, we cannot tell the origins of visible particles: whether the particles 
come from decays of heavier particles or are QCD radiations. Practically, leading two jets 
in pt are used as visible particles to construct Mt2- If they correspond to two "correct" 
particles, namely if each particle is a decay product of each pair-produced particle, Mt2 
behaves as discussed above. However, the leading particles can be decay products of one 
parent particle, and also hard ISR can be one or both of the leading particles. These cases 
are called "combinatorics". 

In many events, Mt2 of the leading particles corresponds to Mt2 of the correct particles. 
For instance, the rate of correspondence is about half for q^^^q^^'^ or ti as shown later. Com- 
binatorics smears the Mt2 distribution. The smearing effect is significant for high Mt2, and 
it is different in each process. 
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Figure 2: Distributions of M'r2 for g'^^^g'-^^ 



gg^(l)^(l) 



— )■ qq^)"-"'^^"' in the left and ti — )■ hhW^W~ 
in the right generated by MadGraph/MadEvent 4.4 [ID] where m^(i) = 912.5 GeV, 
m (1) = SOO.lGeV, and /xq = 211.0GeV. Mj'2 is constructed with the correct two partons 
(green shaded area), the leading two partons (blue solid line), and leading two jets (dotted 
line). Also, events were generated without additional jets, that is, without USR in the 
parton level, and Mt2 was calculated with the correct two partons (red shaded area). The 
distributions are normalized to 1. 



In order to see the effects of USR and combinatorics, we generated the g^^^g*^^^ production 
of the MUED benchmark point and the ti production adding up to one jet in the parton level 
G, and constructed Mt2 with the correct partons and with the leading partons/jets. These 



^The events were generated by MadGraph/MadEvent 4.4 at -^s = 14 TcV. The fragmentation and 
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Mt2 distributions are shown in fig. 2. 

The green (red) shaded area shows Mt2 with (without) an additional jet in the parton 
level using correct two partons: quarks from the direct decay g^^-* — >■ q^'^^^ for q^^\'^^^ and 
h quarks from t — )■ hW for ti. LKPs and neutrinos produce p™*^. For the signal, Mt2 
without USR (red shaded area) is bounded by the mass combination /iq = 210 GeV given 
by f fT3|) . Including USR (green shaded area) Mt2 varies with Pt and exceeds /io as flT^ . and 
theoretically it reachs about 440 GeV with extremely large Pt- For the ti background, Mj'2 
of correct two partons with USR (green shaded rarea) does not exceed rratop as expected. 



Parton level 


qWqW ^ qq^W^W + 0, 1 jet 


ti -^ bbW+W- + 0, 1 jet 


^leadtng ^ ^c^rect 


61.6% 


49.1% 


pleading ^ ^cc^rect 


30.3% 


22.4% 


j^leaatng ^ ^correct 


8.1% 


28.5% 



Table 2: The evaluation of combinatrics for q^^^q^^^ of the benchmark point and ti. Mt2 is 
constructed in the parton level. We compare Mt2 of leading partons and correct partons in 
each event. 

When the leading partons are used for Mt2 (blue solid line), combinatorics smears the 
Mt2 distribution. Mt2 of the leading partons spreads over the endpoint of Mt2 of correct 
partons to roughly double the value at the endpoint as in fig. 2. Table [2] shows the evaluation 
of combinatorics. The leading partons correspond to the correct partons for 60% of the time 
for g^^^g^^-* and for 50% for ti. The smearing effect due to combinatorics is different in 
each process: M^2 *"^ > M^2^^'^^ for three quarters of the combinatorics events of q^^\^^\ 
while M,^2 *"^ < M^2^'^^^ for more than half of the combinatorics events of ti. Therefore, 
combinatorics assists to enhance the signal to background ratio for high Mt2- 

Also, the detector effects are simulated after the fragmentation and hadronization, and 
Mt2 is constructed with the leading two jets (dotted line). The distribution is similar with 
the Mt2 distribution with the leading two partons. 

The dependence on USR makes the signal excess in the high Mt2 region even for the nearly 
degenerate model, and the smearing effect of combinatorics enhances the excess. It can been 
seen that events with high Mt2^ say above 200 GeV, are dominated by the g^^-'g'-^^ signal 
over the ti background. Since for the other SM background processes the parent particle is 
lighter than the top quark, those background events are expeced to have Mt2 lower than 
??7,top. Hence, Mt2 is an effective event selection to search for the nearly degenerate model. 



hadronization were simulated by Pythia 6.4, and the detector effects were simulated by PGS 4 [41]. The 
simulation setup is the same as the simulation described in Sec. 4 except that here only one jet was added 
as the Matrix Element correction. The MLM matching [42] was prescribed. 
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4 Simulation 

4.1 Signal 

Monte Carlo (MC) samples of MUED signal were generated both by a private implementa- 
tion in MadGraph/MadEvent 4.4 (MG/ME) [ID] and an implementation [13] in Pythia 
6.4 [33]. CTEQ5.1L was used as the leading-order (LO) parton distribution function (PDF). 
In the case of MG/ME, the Matrix Element was calculated by HELAS [H], and the frag- 
mentation and hadronization were simulated with Pythia. The effects of jet reconstruction 
and detector smearing were simulated through PGS 4 [H]. 

We consider l/R from 400 GeV to 1.6 TeV in steps of 100 GeV with KR = 10, 20, 30, 
and 40. The remaining parameter, rrih, is set to 120 GeV. The MUED spectrum is simplified 
by neglecting msM for the first KK level. The processes we consider are pair productions 



of the colored first KK particles, g^ ', q^ ' , and Q^ ' . The signal events corresponding to 
luminosities of 5 fb~^ at ^/s = 7 TeV and of more than 10 fb~^ at y/s =14 TeV were 
generated by Pythia. Table [3] shows the production processes which can be generated in 
the Pythia implementation. 



Process 


flavor 


KK gluon + KK gluon 


g + g^ g^^^ + g^^^ 




KK quark + KK gluon 


g + q^g^'^Q'^'^; ^W + ^(i) 




KK quark -|- KK quark 


qr + qj ^ Qf^ + q^^ 


all z,j 
all z,j 


KK quark + KK antiquark 


g + g^QW+QW. qW + qW 
q + q^QW + QW; qW + qW 

q. + q,^Q\'' + qf 

q. + q,^Qf' + Qf^;qr' + qr 

g. + g. ^ Qf^ + QT 


^^3 
i^ 3 



Table 3: Processes of MUED generated by the Pythia implementation [43J. SU{2)l doublet 
and singlet first KK quark are denoted by Q^^-* and g^^-* respectively. 



Since we use the Mt2 dependence on USR, the ISR has an important role. In order to 
reliably evaluate the hard ISR, we considered the Matrix Element correction in MG/ME 
adding up to one jet to the pair productions. The MLM matching [12] was applied to remove 
the overlap between jets from the Matrix Element and ones from the Parton Shower. This 
prescription was demonstrated for the benchmark point of KR = 20 and l/R = 800 GeV. 
The spectrum of this point is listed in table [TJ However it is very time consuming to generate 
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all of the signal MC samples with the Matrix Element correction, so we used Pythia rather 
than MG/ME to generate them for the discovery study. We will show in Sec. 15.21 that 
MC samples generated by ME/ME with the Matrix Element correction have larger excess 
over background than ones generated by Pythia for the benchmark point. Hence, the event 
generation by Pythia is conservative. 



4.2 Background 

MC samples of the SM background, ti, W/Z + jets, Diboson {WW, WZ, and ZZ), etc., were 
produced with MG/ME using the PDF set CTEQ6.1L, and fragmentation and hadronization 
were simulated with Pythia in the same way of the signal. For ti, W/Z + jets, and Diboson, 
up to two partons were added in the Matrix Element and the MLM matching was prescribed. 
The MC samples were detector-simulated through PCS 4. The dominant background pro- 
cesses, ti and W/Z + jets, were normalized to the next-leading-order (NLO) cross section 
consistent with the inclusive dijet analysis of the ATLAS MC study [15] . 

For the sake of comparison with the 41 + EJp'^'''^ analysis, we generated some multilepton 
background processes, such as four leptons through off-shell Z* or 7*. The luminosities of 
generated SM background are more than 2 fb"'^ at ■\/i= 7 and more than 10 fb~^ at i/s= 14 
TeV, respectively. The summary of the background processes is shown in table HI 



Process 


v^ = 7 TeV 


v^ = 14 TeV 


crx efficiency 


crx efficiency 


tt + 0,1,2 jets 


130 pb 


765 pb 


(Vr^/i/) + l,2 jetst 


287 pb 


897 pb 


(Z^/+/-,z/z/) + l,2 jetst 


98.8 pb 


341 pb 


W+W- + 0,1,2 iets 


35.0 pb 


97.5 pb 


W^Z + 0,l,2jets 


15.0 pb 


44.1 pb 


ZZ +0, 1, 2 jets 


4.78 pb 


12.9 pb 


Z*/YZ*/Y -^ 2/+2/- 


29.2 fb 


57.2 fb 


Z + bb 


35.5 pb 


135 pb 


W + bb 


22.5 pb 


55.0 pb 


{Z/Y ^l^l~,vv)+ti 


28.7 fb 


205 fb 


{W -^ Iv) + ti 


37.2 fb 


131 fb 



Table 4: Summary of the SM datasets used in this analysis. -|-n jets are added in Matrix 
Elements with the MLM matching. ''"MC samples of W/Z + jets are preselected as p^^^^^ > 
90 GeV. The cross sections listed are calculated with the preselected MC samples generated 
by MadGraph/MadEvent 4.4. 
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5 Analysis 

5.1 Object selection 

The object selection is that an electron and a muon are required to have pt > 10 GeV 
and |?7| < 2.5 and a jet is required to have pr > 20 GeV and |?7| < 2.5. In order to avoid 
recognizing a shower from an electron as a jet, a jet within AR < 0.2 (Ai? = \/Ar]'^ + A0^) 
from any electron is removed. Charged leptons from hadronic activity also should be removed. 
If an electron and a jet are found within 0.2 < AR < 0.4, the jet is kept and the electron is 
rejected Similarly, if a muon and a jet are found within AR < 0.4, the muon is rejected. 

5.2 Event selection 

First, we require two jets with p^*^* > {100, 20 GeV}, and p™**"* (E™''*) must exceed 100 
GeV. Cuts stronger than these will reduce the MUED signal. At least one charged lepton, an 
electron or a muon, with p^^ > 20 GeV is required in our analysis. The cuts imposed above, 
particularly the requirement of one lepton, are necessary to avoid the QCD background. 
Either {p^^^ > 100 GeV,p7^" > 100 GeV} or p^^ > 20 GeV is used as a trigger. When 
there is only one lepton, we impose a cut of m^^'"^^^" > 100 GeV to reduce the W + jets 
background, where 



^lep,miss _ U U-Pp^iss _ p^^P . ^rni 



Finally, we construct Mt2 with the leading two jets, and impose Mt2 > 200 GeV. To 
summarize our event selection: 

• CUTl: p^' > {100,20 GeV} 

• CUT2: EJf''^' > 100 GeV 

• CUTS: At least one lepton with p^^ > 20 GeV 

• CUT4: If the number of lepton is one, M^^'™'" > 100 GeV 

• CUT5: Mt2 > 200 GeV. 

In order to demonstrate the effectiveness of Mt2 for the MUED, we only use the basic cuts 
1-4 except one on Mt2- The cuts 1-4 are comparable with ones imposed in the ATLAS and 
CMS new physics searches in one lepton + jets + E'^'''^^ with low luminosity at ^/s = 7 TeV 

W7\ . We do not use the Mf,ff cut and the E^p'^^'^'/M^ff cut which are used to extract the 
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signal especially by the ATLAS collaboration in the search for SUSY [13] , where 



jet lepton 



miss 
T 



(16) 



It is common that the A(f)jet,miss cut is applied to reduce events with fake missing due to the 
mismeasurement of jets, but this is not necessary because the later Mt2 cut has a similar 
role ED. 



Process 


CUTl 


CUT2 


CUT3 


CUT4 


CUT5 (Optimal) 


^(i)+^(i) 


MG/ME 


1,028 
937 


832 

757 


119 
108 


62 
63 


25 
22 


PYTHIA 


g^'^ + q'^'VQ''^ 


MG/ME 


9,196 
8,569 


7,218 
6,694 


1,234 
1,344 


675 

731 


241 
223 


PYTHIA 


gW/gw 


MG/ME 


5,315 
4,497 


4,035 
3,276 


863 
690 


508 
436 


148 
84 


PYTHIA 


gW/gw 


MG/ME 


1,444 
1,301 


1,075 
955 


206 
163 


115 
112 


27 
20 


PYTHIA 


Total MUED 


MG/ME 


16,983 


13,160 


2,422 


1,360 


441 


PYTHIA 


15,304 


11,682 


2,305 


1,342 


349 


tt 

W 

z 

W/Z + ti/bb 
Diboson 


426,074 

400,527 

142,368 

1,121 

29,141 


57,533 

97,907 

53,801 

304 

4,482 


23,239 

35,386 

916 

103 

1,335 


5,620 
1,031 

107 
49 

252 


243 
85 
12 
10 
40 


Total Standard Model 


999,231 


214,027 


60,979 


7,059 


390 


Total MUED 
Zb 


MG/ME 


0.05 


0.17 


0.06 


0.78 


4.10 


PYTHIA 


0.05 


0.14 


0.05 


0.77 


3.37 (7.57) 



Table 5: Cut flow for 1 fb"^ at y/s = 14TeV. The MUED benchmark point is {1/R,AR} = 
{800 GeV,20}. The MUED signal generated by MadGraph/MadEvent 4.4 (MG/ME) 
with the MLM matching is normalized to one generated by Pythia. Mt2 > 350 GeV is an 
optimal cut that maximizes the significance Zb = 7.57. 

The cut flow in table [5] shows that the Mt2 cut (CUT5) significantly reduces the SM 
background to a level comparable to the MUED. Since the Matrix Element correction in- 
creases event rates in the high Mt2 region, there remain more signal events generated by 
MG/ME after CUT5 than ones generated by Pythia. Therefore, the signal rate in the fast 
event generation by Pythia is a little underestimated and hence is conservative. 

Fig. [3] shows that the dominant background are ti, W + jets, and Diboson. Background 
events that remain even after CUT5 mainly come from combinatorics. The peak of MUED 
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events is Mt2 < 200 GeV, but the signal events have a long tail which can be understood as 
a result of the variant endpoint due to the wrong test mass discussed in section 13.21 There 
is combinatorics for both the signal and the background, and especially it tends to increase 
Mt2 of the signal events. Combinatorics help to enhance the signal excess. As a result, we 
can successfully extract the signal from the background even based on jets. 



M-j-2 distribution of background at 14TeV 




SM background 

tt 
W 
Z 

Diboson 

W/Z+bb/tt 



-i^H 



LJh 



Af J, distribution of MUED signal at 14TeV 

I/R=700GeV, AR=20 
I/R=800GeV, AR=20 
I/R=900GeV, AR=20 
SM background 




100 200 300 400 500 600 700 800 900 1000 



100 200 300 400 500 600 700 800 900 1000 



Figure 3: Distributions of Mt2 of the leading two jets after CUT4 for each SM background 
in the left, and total SM background and MUED signal points of l/R = 700, 800, 900 GeV 
and AR = 20 generated by Pythia in the right. 



5.3 Discovery potential 

For the study of the discovery potential it is necessary to take systematic uncertainties into 
account in addition to statistical uncertainties. We use the significance Zb |1H], which is 
provided by the ROOT library [19], using the same approach as in the ATLAS discovery 
study of the SUSY |15]. Zb is calculated using a convolution of a Poisson and a Gaussian 
term to account for systematic errors. For the backgrounds except those from QCD, a 
reasonable estimate of the systematic uncertainty is ±20%. The discovery potential is studied 
by finding the optimal Mt2 cut (in step of 50GeV) to maximize the significance Zb- We define 
"discovery" when Zb > 5 and more than 10 signal events remain after the cuts. 

In order to compare the Mt2 analysis in the multijet + lepton mode with the previously 
studied 4/ + i?™** analysis, we also check the discovery potential in 4/ + i?™** using the 
same MC samples and using the same definition of discovery. In the 41 + E'^*^'* analysis, the 
following cuts are imposed [Il]:(l) four isolated leptons with p^^ > {35,20, 15, 10 GeV} are 
required, (2) E"™*^ > 50 GeV, and (3) an invariant mass Mn for all possible pairs of opposite 
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sign same flavor leptons and remove events if \Mii — mz\ < 10 GeV to reduce background 
from the Z boson. Tlie estimated background from our MC samples is 10 events/100 fb~^. 
The fake leptons should be considered to evaluate the background level of 4/ + EJp'''^^ more 
appropriately, but the fake leptons are not considered since they are not important for our 
analysis based on jets. 



> 
O 



1400 
1300 
1200 
1100 
1000 

900 ^ 

800 

700 

600 

500 

400 
10 



LHC 1fb-' limit 




^^ Mt2 lOfb' 14TeV 

41epton lOfb' 14TeV 

—^- Mj^ Ifb' 14TeV 
41epton Ifb' 14TeV 



M^j 2fb"' 7TeV 



Figure 4: Discovery potential of the MUED with 1 fb~^ and 10 fb~^ at ^/s = 14 TeV in the 
4/ + EJp^'^'^ analysis and the M7-2 analysis and discovery potential of 2 fb~^ at ^/s = 7 TeV 
only in the Mt2 analysis. For a given luminosity, the parameter region below the line will be 
discovered. 



The spectrum is more degenerate for smaller AR, which is more difficult for discovery 
in general. Note that for fixed 1/R, the MUED with smaller AR has a larger cross section 
simply because the KK gluon and the KK quark become lighter as in Eqs. ([7]) and (|8]). Fig. 4 
shows that the discovery potential does not vary with changing AR in the Mt2 analysis. The 
first run of LHC at ^/s = 7 TeV will have an integrated luminosity of 2 fb~^ and so it can 
discover up to 1/R ~ 500 GeV. However this parameter region was already excluded by the 
ATLAS multijet+£'™'^* analysis with 1 fb~^, as mentioned in Sec. 12. 4[ The second run at 14 
TeV will discover up to l/i? ~ 1 TeV with 1 fb^^ and 1/i? ~ 1.2 TeV with 10 fb'^. 

In the 4/ + E^''"^ analysis, the discovery reach at 14 TeV is 1/i? ~ 700 GeV with 1 fb"^ 
and 1/R ~ 1.2 TeV with 10 fb~^ for 20 < AR < 40, but the sensitivity is very low for 
Ai? = 10 : the discovery reach is only 1/R = 400 GeV with 1 fb"^ and 1/R = 800 GeV with 
10 fb-i. 
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The result shows that our Mt2 analysis improves the discovery potential. In particular, 
the improvement is so significant for the most degenerate parameter KR = 10 that the 
discovery potential improves from 1/R = 400 GeV to 900 GeV. 

6 Discussion and Conclusion 

We have improved the discovery potential of the MUED. We include the background sys- 
tematic uncertainties. The previously studied multilepton channels, such as 4/ + E'™**, are 
clean channels but only sensitive to the limited number of processes. On the other hand, 
we utilize the multijet + lepton channel, and it is accessible to the most processes of the 
MUED. Although this channel is statistically advantageous, it difficult to extract the MUED 
signal events from backgrounds. However, we succeeded in extracting the signal from the 
background by applying Mt2 for the event selection. This is because the signal events can 
have large values of Mt2 thanks for the dependence on USR with the wrong test mass, and 
because the combinatorics effect enhances the signal excess. 

The simulation result indicates that the discovery potential is significantly improved by 
using Mt2 cut. The most important achievement is that our Mt2 analysis has a much greater 
sensitivity for the MUED for the most degenerate mass spectrum than the 4/ + £'™** analysis. 
For 1 fb~^ at ^/s = 14 TeV, the discovery potential increases by 500 GeV in terms of 1/R 
for the fixed AR = 10. Therefore, the Mt2 analysis is particularly effective for the case of a 
highly degenerate spectrum. 

The relic abundance of the Dark Matter implies that the very high value of 1/R ~ 1.5 TeV 
is favored, while our analysis reaches only up to 1/R ~ 1.2 TeV with 10 fb~^ at ^/s = 14 TeV. 
However, there is possibility to improve the discovery potential further. We considered the 
Matrix Element correction and the NLO correction to the background. If those corrections to 
the signal are included, the Matrix Element correction enhances the event rate for high Mt2 as 



shown in Sec. 15. 2[ and the NLO correction increases the cross section. Also, additional cuts, 
such as the 6-jet veto and E^^'^'^/Mefj, could enlarge the discovery potential of the MUED, 
although we did not introduce them in order to emphasize the effectiveness of Mt2- The 
fe-jet veto must be particularly effective because it significantly reduces the ti background. 
Including these effects the discovery reach will be improved, and the Dark Matter favored 
MUED may be searched with O(IOO) fb"^ 

The search using Mt2 has a model independent aspect. The benefit of the Mj'2 analysis is 
significant when there is large USR. We might expect large USR because a common feature 
of many new physics models is production of heavy colored particles associated with hard 
ISR contained in USR. This analysis is also applicable to other models: SUSY with R parity 
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or Little Higgs model with T parity with a nearly degenerate mass spectrum. 

For our analysis, we required one lepton to avoid the QCD background, but this is not 
always necessary. If the well-understood QCD MC samples are prepared, we can study 
the discovery potential based on the Mt2 analysis in the inclusive dijet channel, which is 
statistically advantageous and is accessible to models that rarely emit leptons. 

Because of the significant improvement of the search for MUED by using Mt2i this work 
implies that the LHC will be able to have better sensitivities to nearly degenerate models. A 
more general analysis with a nearly degenerate spectrum is left to a future work. 
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